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Abstract 

We present the results of molecular dynamics simulations of the Gay-Berne 
model of liquid crystals, supercooled from the nematic phase. We find a glass 
transition to a metastable phase with nematic order and frozen translational 
and orientational degrees of freedom. For fast quench rates the local structure 
is nematic-like, while for slower quench rates smectic order is present as well. 
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The study of nematic glasses [|| is richer than the corresponding study of isotropic glasses 
due to the presence of orientational degrees of freedom in the former systems. Assuming 
that we supercool a liquid crystal starting from its nematic phase (rather than the isotropic 
phase) we expect a nematic glass to have long-range orientational order like an equilibrated 
nematic as well as frozen density and director fluctuations. Orientationally this glassy phase 
is similar to a mixed magnetic phase where both ferromagnetic and vector spin glass order 
coexist 0. However, unlike a spin glass the nematic glass is a nonequilibrium phase of matter 
and the molecular translational degrees of freedom freeze as well at the glass transition. 

In this Letter we study the formation of a nematic glass using molecular dynamics (MD). 
We model the liquid crystal using the Gay-Berne (GB) potential 0] which is an anisotropic 
Lennard- Jones potential. Previous molecular dynamics studies have indicated that the Gay- 
Berne potential exhibits a rich phase diagram ||^-|^ and the principal dynamical features 
HTHTOj of real liquid crystals. The GB potential is given by, 



where ui,U2 are unit vectors giving the orientations of the two molecules separated by the 
position vector r. The parameters e(ui,U2,r) and cr(ui,U2,r) are orientation dependent 
and give the well depth and the intermolecular separation where U = respectively. The 
well depth is written as. 
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The range parameter a{ui, U2, r) is given by 

a(ui,U2,r) = ao<l- -< 



2 I 1 + x(ui ■ U2) 

2 s ^ -1/2 



[r ■ ui - r ■ U2 



1 - x(ui • U2) 
The shape anisotropy parameter x is given by 



(5) 



X = {(a>,)2-l}/{(a,M)2 + l} (6) 

where cxe and as are the separation of end-to-end and side-by-side molecules respectively. 
The parameter x' is given by 

x' = {l-{ee/es)'/'}/{l + {eJes)'/'} (7) 

The ratio of the well depths for end-to-end and side-by-side configurations is Sej^s- 

We investigated glass formation in the Gay-Berne fluid using a constant-pressure, 



constant-temperature MD method |Tl|], which allows the volume to change as a function 
of the temperature and pressure of the system. The edges of the cell were allowed to vary 
independently, but the orthogonal shape was maintained. We simulated a system of = 864 
particles. We chose Oejos = 3, Ce/cs = 5, and u = 1 and = 2 as in the original work of Gay 



and Berne. The moment of inertia was chosen to be Amal, as in ref. |T^. We used periodic 
boundary conditions, and cut off and smoothed the potential at 3.8o"o- The equations of 
motion were solved using the leap-frog algorithm with an integration time-step At* = 0.001 
in dimensionless units (At* = At(mcrQ/eo)^^^^, where m is the mass of a molecule ). The 
initial configuration was chosen in the nematic phase at the dimensionless temperature 
T*(= Aj^T/eo) = 1.2 and pressure P*{= PaQ/eo) = 5.8 ||12|. The nematic order parameter 
at this point was 5* = 0.75. The temperature was then reduced in a sequential fashion 
allowing the system to relax at each temperature for a number of time steps depending on 
the quench rate. The quench rate is defined by the ratio of the change in dimensionless 
temperature to the number of timesteps between two succeeding temperature reductions. 
Our discussion will focus primarily on simulations where the temperature was lowered in 



decrements of AT* = 0.1 every 100 iterations, corresponding in real units to a quench rate 
approximately equal to 10^^ K/s. Rates faster and slower than this one will be discussed 
at the end. After the system was cooled to its final temperature, it was allowed to anneal 
and various structural and dynamical properties were measured. Fig. |^ shows the potential 
energy as a function of the annealing time at four different final temperatures. The breaks in 
the curves corresponding to temperatures T* = 0.4 and T* = 0.5 suggest that crystallization 
may have occurred. Further proof of crystallization can be found by examining pair distribu- 
tion functions after the breaks take place. In Fig. ^ we show parallel, perpendicular (relative 
to the director) and orientationally-averaged pair distribution functions at T* = 0.4 after 
an annealing run of 120,000 iterations. From the plot of the parallel distribution function 
we see that a smectic density wave has been formed parallel to the director. The peaks of 
the perpendicular distribution function show the hexagonal crystalline ordering within the 
layers. The splitting of the second peak of the orientationally-averaged distribution func- 
tion is typical of a crystalline phase. The nonzero value of the perpendicular distribution 
function (which is averaged over the layers) at r = shows the strong correlation between 
the layers. At T* = 0.3 and below, a different picture is observed even after annealing. 
The nematic order parameter reaches a plateau at long times (see Fig. with average 
values lower than the corresponding values in the equilibrium phase at the same tempera- 
ture {S = 0.878, 0.920 in the quenched system, S = 0.964, 0.962 in the equilibrium system 
for T* = 0.2,0.3 respectively). The pair distribution functions for a system quenched to 
T* = 0.2 and then annealed for 120,000 iterations are shown in Fig. ^ There is no smectic 
layer formation and no translational order perpendicular to the director, i.e. structurally 
the quenched system is nematic-like. 

We now show that the system is frozen in this nematic state both translationally and ori- 
entationally. The dynamical properties of the system are characterized by the translational 
diffusion coefficients and orientational relaxation times. Diffusion coefficients parallel and 
perpendicular to the director are evaluated from mean-squared displacements as follows: 
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D|| = lim - < (r|| (t + to) - ry (to))' > (8) 
= lim ^ < (r^(t + to) - rx(to))2 > . (9) 

t^oo t)t 

Both diffusion constants in the glassy phase are very small, about 100-1000 times smaller 
than corresponding values in the nematic phase which are of order 0.1 in dimensionless units. 
However, the mean-squared displacement remains a monotonic function of time, unlike the 
crystalline phase where the mean-squared displacement oscillates about the equilibrium 
positions. The ratio of the two diffusion constants is approximately 4 in the glass phase, 
comparable to values in the nematic phase, and very different from the equilibrium values in 
the smectic and crystalline phases where the ratio is less than one and approximately equal 
to one respectively. 

To show that the system freezes orientationally we use the spin-glass analogy and com- 
pute an Edwards-Anderson type |]14| of correlation function: 

1 ^ 

Cit) = -T- E{< + ^o)u.(to) > - < Ui{t + to) >< U,(to) >}. (10) 

i=l 

The orientational relaxation time is given by the slope of the logarithmic plot of C(t) versus 
time t and is approximately 2 orders of magnitude larger in the glass phase compared with 
the nematic phase, as shown in Table |. Thus, our evidence suggests that a glass phase 
has formed at this temperature with simultaneous freezing of translational and orientational 
degrees of freedom. 

Finally we consider the effects of different quench rates. If we quench the system instan- 
taneously from T* = 1.2 to T* = 0.2 we obtain a glass. However, if the glass is then annealed 
for several thousand time steps crystallization occurs. This shows that the glass becomes 
more unstable with increasing cooling rate, as seen in simulations of isotropic Lennard- Jones 
systems [|T^. On the other hand, with a "slow" cooling rate corresponding to a tempera- 
ture reduction of AT* = 0.1 every 1000 iterations we find a stable glassy phase. However, 
the slower quench yields a different structure than the one shown in Fig. § because of the 
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smectic phase that intervenes between the nematic and crystalhne phases under equihbrium 
conditions. At the slower coohng rate the system starts to form layers typical of the smectic 
phase, but does not crystallize. There are no correlations between the different layers and 
order within the layers is not perfect. These effects can be seen in Fig. ^ where we show the 
positions of the molecules in three adjacent layers of the MD cell when viewed from above 
for the "slowly" quenched system at T* = 0.2. Fig. ^ shows a similar view in a system which 
is allowed to reach thermal equilibrium at the same temperature and a crystalline phase is 
formed. 
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FIGURES 

FIG. 1. The potential energy in dimensionless units as a function of annealing time (also in 
dimensionless units) for a Gay-Berne fluid quenched from T* = 1.2 to the temperatures indicated. 
The quench rate corresponded to a reduction in temperature of 0.1 every 100 iterations. Note the 
absence of a break in the curves corresponding to the temperatures T* = 0.2 and 0.3, unlike the 
two higher temperatures. 

FIG. 2. The pair distribution function versus distance in dimensionless units (r* = r/ao) 
after a quench to T* = 0.4 with the same quench rate as in Fig. |l[ The distribution functions 
were evaluated after an annealing run of 120,000 iterations, and averaged over the subsequent 
120,000 iterations. The solid curve is the orientationally averaged function. The dashed curve 
is the distribution function parallel to the director (averaged over the perpendicular direction) 
plotted against distance parallel to the director. The dashed-dotted curve is the distribution 
function perpendicular to the director (averaged over the parallel direction) plotted against distance 
perpendicular to the director. 

FIG. 3. The nematic order parameter as a function of annealing time for quenches to the 
temperatures indicated. The quench rate is the same as in Fig. 0. Note the lower value of the 
order parameter for the lowest of four temperatures shown. 

FIG. 4. The pair distribution function versus distance in dimensionless units after a quench to 
r* = 0.2 with the same quench rate as in the previous figures. The distribution functions were 
evaluated after an annealing time of 120,000 iterations, and averaged over the subsequent 120,000 
iterations. The line style of the three curves is as in Fig. ^. Note the absence of smectic order in 
the parallel distribution function (the lowest curve), and the absence of two well- formed subsidiary 
peaks in the orientationally-averaged distribution function (the solid curve). Compare with Fig. 
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FIG. 5. The locations of the molecules in the three innermost smectic layers of the MD cell 
(there are fifteen layers in total) for a system quenched to T* = 0.2 at the "slow" cooling rate 
and then annealed for 200,000 iterations. The triangles correspond to molecules in the middle 
layer of the three layers shown. The boxes and circles correspond to the two layers which flank 
the middle one. Note the absence of correlations in molecular positions between layers, and the 
imperfect hexagonal order within each layer. Compare with Fig. |6| which shows a system at the 
same temperature in thermal equilibrium. 

FIG. 6. The locations of the molecules in the three innermost smectic layers of the MD cell for 
a system in thermal equilibrium at T* = 0.2. The notation used is the same as in Fig. ^. In this 
run, which produced an equilibrated crystalline phase, the potential was truncated perpendicular 
to the director at a distance of 2.5. Two side- by-side molecules at that separation have the same 
potential energy as two end-to-end molecules separated by the spherical cutoff distance of 3.8. 
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TABLES 

TABLE I. The orientational relaxation time constant r and the nematic order parameter S at 

several representative temperatures for a nematic in thermal equilibrium and a glass quenched at 
two different rates. The fast quench rate corresponds to a temperature reduction of 0.1 every 100 
iterations and the slow rate to the same reduction every 1000 iterations. The glass formed at the 
faster quench rate has less nematic order and significantly longer relaxation times. 







r 


S 


Nematic, T* = 1.0 




5.1 ±0.2 


0.84 ±0.01 


Nematic, T* = 1.1 




5.2 ±0.2 


0.79 ± 0.01 


Glass, slow quench, T* = 


0.2 


394 ±6 


0.888 ± 0.004 


Glass, slow quench, T* = 


0.3 


152 ± 1 


0.904 ± 0.007 


Glass, slow quench, T* = 


0.4 


48 ±2 


0.969 ± 0.002 


Glass, fast quench, T* = 


0.2 


676 ±9 


0.878 ± 0.004 


Glass, fast quench, T* = 


0.3 


93 ±1 


0.920 ± 0.004 


Glass, fast quench, T* = 


0.4 


100 ±2 


0.961 ± 0.002 
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